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We study effect of inhomogeneities of the lattice in the 
S — 1/2 bond-alternating chain by using a quantum Monte 
Carlo method and an exact diagonalization method. We 
adopt a defect in the alternating order as the inhomogeneity 
and we call it "bond impurity". Local magnetic structures 
induced by the bond impurities are investigated both in the 
ground state and at very low temperatures. The local mag- 
netic structure can be looked on as an effective 5 = 1/2 spin 
and the weakness of the interaction between the local struc- 
tures causes the quasi-degenerate states in the low energy. We 
also investigate the force acting between bond impurities and 
find that the force is generally attractive. We also study the 
dynamical property of the local magnetic structure. While 
the local magnetic structure behaves as an isolated 5 = 1/2 
spin in the response to a time-dependent uniform field, it is 
found to be robust against the effect of a random noise applied 
at each site individually in the sweeping filed. 



I. INTRODUCTION 

In the low dimensional quantum spin systems, it has 
been turned out that the quantum mechanical effect, in 
particular the mechanism of the singlet pair formation, 
plays an important role for the ground state spin config- 
uration . There have been found a wide variety of 
peculiar arrangements of the singlet pairs in the ground 
state which brings new types of singlet (or nonmagnetic) 
ground state phases such as VBS ||, RVB [|]||, etc. As 
general characteristics of these systems, a finite energy 
gap exists between the ground state and excitation states 
and the correlation length in the ground state is finite. 
If a lattice has some inhomogeneities such as impurity 



sites, defects of periodicity, and edge points, etc., then lo- 
cal magnetic structures are induced around them ||-^). 
Property of such induced moments is one of the most 
interesting current topics. 

In this paper we will study the interaction between 
the induced moments and also dynamical properties of 
such moments in the bond-alternating Heisenberg anti- 
ferromagnetic (HAF) chain. This model is one of the 
simplest systems of the singlet ground state which con- 
sists of singlet states at strong bonds. The magnetic 
susceptibility and specific-heat of this model was stud- 
ied in detail as a function of the ratio of the alternation 
coupling constants by Duffy and Barr |u|. Various ma- 
terials corresponding to this model have been studied ex- 
perimentally, i.e., aromatic free-radical compounds [ pd| , 
Cu(N0 3 ) 2 -2.5H 2 ||, and (VO) 2 P 2 7 @, etc. The 
ratio of the alternation coupling constants for theses ma- 
terials has been estimated comparing with the theoretical 
results 11 1 . This system has been also extensively stud- 
ied concerning with the spin-Peirels transitions |l~4] |l(| . 

Here we consider a defect of the alternating or- 
der of the bonds, such as • ■ ABABAABAB- • • or 
• • -ABABBABAB- • •, where A or B denotes the strength 
of bonds (the magnitude of exchange constant) . At this 
defect the singlet dimer state cannot be formed there. 
First we study how a magnetic structure appears at these 
positions. In this paper we call such a defect of the alter- 
nating order "bond impurity". Furthermore, the induced 
magnetic structure depends on whether the edge bond is 
A or B. Bond impurity effects in the uniform HAF have 
been investigated in our previous paper 17 and we com- 
pare the role of the bond impurity in the present model 
with them. 

Next, we focus on what kind of interactions act be- 
tween these locally induced magnetic structures. We 
study the force by investigating the dependence of the 
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ground state energy on the distance between the impu- 
rities for various types of configurations of defects. The 
force between the bond impurities has been also studied 
for the uniform HAF There we found that the at- 

tractive force acts. In the present model we again find 
only attractive force regardless of the types of configura- 
tion. 

Due to the finite correlation length of the present 
model, the induced moments behave almost indepen- 
dently. Thus we can regard such a local moment as 
a microscopic moment. Because recently the technol- 
ogy in microscopic processing makes remarkable progress 
and the analyses in microscale or nanoscale phenomena 
have become possible, the quantum phenomena in mi- 
croscale or nanoscale have received much attention. For 
example, the resonant tunneling phenomena have been 
observed in recent experiments on high-spin molecules 
(Mni2) Jl8|-p0t and the concept of quantum tunneling of 
the magnetization (QTM) has become a topic of interest. 
Recently the adiabatic motion of spin 1/2 was observed 
in V15 and the effects of thermal disturbance on the dy- 
namics has been discussed pl[ . We have studied the 
tunneling dynamics from viewpoint of the nonadiabatic 
transition 

Under a time-dependent magnetic field, how does that 
local magnetic structure behaveH We study the response 
of the magnetization to a sweeping magnetic field. Al- 
though the induced moment consists of several spins, 
the total magnetization behaves as a single spin when 
a uniform field is applied. Under a sweeping magnetic 
field the behavior of magnetization is described by the 
Landau-Zener-Stiickelberg (LZS) formula |p6|~p8||. In re- 
alistic situations, noise disturbs the simple behavior of 
LZS formula |p9|j30f . Furthermore, if the noise acts inde- 
pendently at each site, the behavior of induced moment 
(a cluster of spins) is expected to be different from that 
of isolated single spin. We investigate effects of such in- 
dividual noise at each site on the dynamical properties. 

This paper is organized as follows. In the next section, 
we explain briefly the method used in this study. In Sect. 
[II, effects of bond impurities on the magnetic structure 
in bond-alternating chains are studied. In Sect. |v|, we 
study the force between bond impurities. In Sect. [v|, 
we investigate the response of the magnetization to a 
sweeping field. Sect. VI is devoted to the summary and 
discussion. 



II. MODEL AND METHOD 



where J^i+i changes alternately among Ji (a strong 
bond) and J 2 (a weak bond). Here we consider defects of 
the alternation which cause inhomogeneities of the lat- 
tice (bond impurities). To study low temperature prop- 
erties of the model, we mainly use the loop algorithm 
of the continuous time quantum Monte Carlo method 
(LCQMC) with a method of specification of the 

magnetization M z This method overcomes the 

problem of long autocorrelation in Monte Carlo update 
and allows us to study systems at very low temperatures. 
In the present work, we performed 10 5 Monte Carlo steps 
(MCS) for getting equilibrium of the system and 10 6 
MCS to obtain quantities in the equilibrium state. Here 
a MCS means an update of the whole spins. 

In order to study the dynamical response of spins un- 
der a time-dependent external field, we investigate the 
following system: 

n = Y J J i,i+i s i ■ s *+i + 2r^5f - 5^(2T(f) + h(t))S*, 

i i i 

(2) 

where Hit) is a time-dependent external field and hi(t) 
is a random noise applied to each site individually. T 
is the transverse field which represents terms for quan- 
tum fluctuation of M z . The time evolution of the state 
is obtained by the time-dependent Schrodinger equation 
(TDSE) 



*n^|*(*)> =«!*(*)>> 



(3) 



where \ ^(t)) denotes the wave function of the spin system 
at time t. We set h = 1. Equation (|3|) is solved using the 
fourth-order fractal decomposition ]33[ , 



-un 



[S 2 {-itp 2 )] 2 S 2 (-it{l - 4p 2 ))[S 2 (-itp 2 )} 2 , (4) 



where S 2 (x) = e ^n 1 /2 e xH 2e xHi/2 with p2 = (4 _ 43)-!, 
putting 

T~ti = J]] Ji,i+l&i • Sj+i + 2r , 

i i 
i 

As the initial state, we set the applied field to its min- 
imum value H(t = 0) = — Hq < 0, and put the system 
to be the ground state for this field. Then we sweep the 
filed as 



The Hamiltonian treated in the present paper is given 

by 

T~t = / ] Jj,i+iSj • Sj+i, (1) 

i 

where S, ; = (Sf,Sf,Sf) are the S — 1/2 spin opera- 
tors at the site i. We study the bond-alternating chain 



H(t) = -H + ct. (5) 

Beside this sweeping field, we provide a random noise 
with an exponential-decaying autocorrelation function 
{hi(t)} by a Langevin equation (Ornstein-Ulenbeck pro- 
cess), 

h(t) = - 7 h(t) + r,(t). (6) 
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Here rj(t) is a white gaussian noise, 

(r,(t)) = and (v(O)v(t)) = A 2 S(t), 



where A is amplitude, 7 is damping factor. Thus ob- 
tained random process hit) has the properties 

(h(t))=0 and (h(0)h(t)) = ^-exp(--), (8) 

Z7 T 

where r = 1/7. The dynamics of magnetization is ob- 
tained as 



M(t) = <*(*)! £W*)>- 



(9) 



If the magnetic field changes rather slowly, i.e., the sweep 
rate is rather small, ^(t)) changes adiabatically. In this 
case the state stays in the ground state of the system 
for the current field H(t), and the magnetization follows 
the value of the ground state. When the sweeping rate 
becomes large the system cannot follow the change of the 
field completely, and the nonadiabatic transition occurs. 
The probability to stay in the ground state was given by 
LZS formula as 



p = 1 - exp( ). 



(10) 



When the effect of noise does not become negligible, the 
noise disturbs the quantum process and this probability 
changes [p9[| . 



III. BOND IMPURITY IN BOND-ALTERNATING 
CHAINS 



Here we take the strong bond to be Ji — 1.3 and the 
weak bond to be J 2 = 0.7. For this set of bonds, the cor- 
(7) relation length is estimated as £ ~ 0.82 (see Appendix). 
Here we take &b as a unit of energy (fee = 1)- These four 
models have odd number of spins and their ground state 
is doublet according to the Lieb-Mattis theorem 34 1. We 
performed simulations at T = 0.01 in the M z — 1/2 space 
to study magnetic structures in the low energy state. The 
magnetization profiles {rrii} of (a)-(d) are drawn in Fig. 
[l], where m, L = (S-) and ( ) denotes the canonical aver- 
age at a given temperature. A magnetization is induced 
locally around the impurity. First we consider the cases 
where the bonds at edges are strong, i.e., the model (a) 
and (d) . If we allocate a singlet pair at each strong bond, 
neighboring two strong bonds remain at the center in the 
model (a), while one site remains in the model (d). The 
magnetization of M z = 1/2 is assigned in the remain- 
ing part at the center of the lattices to induce a local 
magnetic structure. Because the edge bonds are strong, 
no magnetization is induced at the edges. Figures 0(a) 
and (d) are considered to describe well the magnetization 
profiles of the ground state since these are gapful systems 
without quasi-degenerate states within a subspace with 
a fixed magnetization (i.e., M z = 1/2). We find that M z 
is distributed only into the ±1/2 space in the simulations 
at this temperature. 

In the models (b) and (c), when we allocate singlet 
state at the strong bond, there are three positions for 
magnetic moments, i.e., the center and both edges. In- 
deed in the both models (b) and (c), magnetizations are 
induced around the impurity and both edges as shown 
in Fig. 1. In order to study the distribution of magneti- 
zation in the lattice, we introduce the summation of the 
magnetization per site from the left edge site 



We investigate magnetic structures in the system (|l|) 
with a bond alternation ■ • • Ji J\ J-x • • •, where Ji > J 2 . 
We study effects of a defect in the alternation, such as 
• ■ • J 1 J 2 Ji J 2 J\ Ji J 2 J 1 J 2 J\ ■ ■ ■■ Beside this defect, mag- 
netic structures can be induced at the edges as naturally 
understood from the VBS picture, i.e., an unpaired spin 
causes an induced magnetization (see the figures). Thus 
we study the following four systems of the bond configu- 
rations: 

(a) a chain of 63 sites where the two strong bonds are at 
the center and both edges terminate with a strong bond, 
( Ji J2 • • • J1J2 J1J1 J2J1 ■ • • J2J1) 

(b) a chain of 65 sites where the two strong bonds are at 
the center and both edges terminate with a weak bond, 
(J2 J\ ■ ■ ■ J1J2 J1J1 J2J1 ■ ■ ■ J1J2) 

(c) a chain of 63 sites where the two weak bonds are at 
the center and both edges terminate with a weak bond, 
(J2 J\ ■ ■ ■ J2J1 J2J2 J1J2 ■ ■ ■ J1J2) 

and 

(d) a chain of 65 sites where the two weak bonds are at 
the center and the edges terminate with a strong bond 
(Ji J2 • • • J2J1 J2J2 J1J2 ■ ■ ■ 'hJi)- 



(11) 



We show this quantity for the model (b) in Fig. ||. There 
the values of the left plateau and right plateau are 0.166 
and 0.333, respectively. From this figure we find a spin 
1/6 locates at each local structure. This deceptive frac- 
tional magnetization is considered to come from mixing 
of states. 

Looking on the local magnetic structure as an effective 
5 = 1/2 spin interacting by an effective exchange </, this 
system is modeled by a three-site Heisenberg model Ti. = 
JSi • S2 + JS 2 • S3 . In the M z — 1/2 space the eigenvalues 
are 



E x = - J, E 2 = 0, and E 3 = J/2. 



(12) 



The corresponding eigenvectors are denoted by \<f>i), (i = 
l,2,and3). The expectation values of magnetization of 
spins in each state are 

(fclfljlfc) = (ik\SS\<h) = 1/3 and (^ilSfl^) - -1/6, 
(fclSflfc) = (hlSllh) = and (MSIlfo) - 1/2, (13) 
{<h\St\h) = (fo\SI\fo) = <&|Sf |03> = 1/6. 
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In the gapped spin system the correlation function de- 
cays exponentially and the effective coupling between the 
induced moments is expected to be very small i.e., 
J -C 1. Thus these states are considered to be almost de- 
generate even at this temperature (T = 0.01), although 
there are energy gaps of order O(J) between the state 
\fa), \<h), an d 1^3)- I n sucn a case the three states ap- 
pear in the equal probability, and the expectation value 
of magnetization is given by an equal-weight average in 
the three state. That is, (Sf), and (S§) are given 

by (1/3 + + l/6)/3 = 1/6, (-1/6 + 1/2 + l/6)/3 = 1/6, 
and (1/3 + + l/6)/3 = 1/6, respectively. These values 
correspond to the observed deceptive fractional magneti- 
zation. 

In order to confirm the above modeling we perform the 
following two investigations. First, we investigate a short 
chain of L = 21 by an exact diagonalization method in 
order to check that the ground state of the type (b) is rep- 
resented by \fa) of the three-spin model. Here we choose 
a chain of a shorter correlation length because the length 
of the chain is short. Namely, we set J\ = 2 and J 2 = 0.5, 
where £ -C 1 . Fig. ^ shows the summation of the magne- 
tization from the left edge site (Eq. (jlT|)) for this model. 
The local magnetization around the left edge is about 
1/3 which corresponds to {fa\S\\fa). This causes the 
left plateau. The local magnetization around the impu- 
rity at the center is about —1/6 ((fa\S%\fa)). Therefore 
the right plateau is 1/6 ((fa \Sf\fa) + (fa\S$\fa)). Fi- 
nally adding the local magnetization around the right 
edge (1/3), Ai z (21) terminates at 1/2. Thus the ground 
state \fa) represents well that of the type (b) model. 

Second, to confirm the quasi-degeneracy in the model 
(b) , we check the distribution of M z in the Monte Carlo 
simulation, which is shown in Fig. 0. Here the distribu- 
tion for M z = 3/2 is about 12.5 %. Because there are 
three states in the M z = 1/2 space and one state in the 
M z =3/2 space in the three-spin model, this distribu- 
tion indicates that these four states are equally populated 
and that T = 0.01 is much higher than the energy gaps 
between these four states. In the three-spin model the 
eigenvalue of the state of M z — 3/2 is J/2, while the 
eigenvalues of the state of M z = 1/2 are — J, 0, and J/2 
(Eq. @). In the M z = 3/2 space we observe S =1/2 
moment at each local structure which represents the state 
I +++) in the effective three-site model. The values {rrii} 
are almost three times as large as those of Fig. [j] (b). In 
principle we can obtain the energy gap from the temper- 
ature dependence of the distribution B. However it is 
too small to detect here. 

We find a similar scenario for the model (c). Thus we 
conclude that in bond-alternating systems, local mag- 
netic structures are induced by a bond impurity or weak 
edge bonds as an effective S — 1/2 spin and they behave 
almost independently. 

Now let us examine more detailed structures of the lo- 
cal magnetic structures. In the case (a) a negative mag- 
netization appears at the middle site, while in the case 



(d) a positive magnetization appears there. The interac- 
tion of the three spins at the center of the model (a) is 
approximately represented by the three-spin model with 
the strong bonds. Here the magnetizations at the center 
of the model (a) are distributed as about (1/3, —1/6,1/3). 
On the other hand in the model (d) a spin at the center 
is isolated from the others. 

The local magnetic structures in the present model are 
well isolated. Therefore we can locate such magnetic 
structures as we desire. In Fig. ^ we show a magnetiza- 
tion profile which has five positions of induced structures 
in a chain of L = 101. This system has three J 2 J 2 defects 
and weak edge bonds. This configuration is obtained at 
a very low temperature (T = 0.01) in the space of the 
total magnetization M z = 1/2. Each moment almost 
behaves independently. Thus magnetization at each lo- 
cation is given by an average over all the possible states. 
In the present case, each local moment has a deceptive 
magnetization 1/10. Furthermore we confirmed that the 
distribution of magnetization P(M Z ) is given by a bino- 
mial distribution 



25 (5-2|AG|)!(2|M,|)!" 



IV. FORCE BETWEEN TWO DEFECTS IN AN 
ALTERNATING CHAIN 

When a bond-alternating system has impurities, what 
kind of interaction does exist between them, attractive 
or repulsive? We study the force between the bond im- 
purities in this section. We estimate the ground state 
energies of the spin system in various fixed configura- 
tions of bonds and compare the ground state energies as 
a function of the distance between the bond impurities. 

In the alternating chain (• • • J\ J2 Ji J2 • • • ), the system 
may have a pair of defects by shifting a position of a 
strong bond by one. 

■■■ JiJ 2 JiJiJ 2 J 2 JiJ2---. (15) 

If we shift the position furthermore, the system has a 
configuration 

■ ■ ■ J1J1J2J1J2J2 ■ ■ ■ , (16) 

etc. We study dependence of the energy on the distance 
(A a ) between the positions of Ji Ji and J 2 J 2 . We define 
A a = for the case of no defect, A a = 1 for the con- 
figuration jl5|), A a = 2 for the configuration ( |l6| ) and 
so on. In Fig. |^ (a) we plot the ground state energy 
as a function of A a obtained by an exact diagonaliza- 
tion for L = 24 with Ji = 2 and J 2 = 1 in the periodic 
boundary condition (PBC). The ground state energy be- 
comes larger as A a becomes larger. Therefore we find an 
attractive force between the impurities (Ji Ji and J 2 J 2 ). 
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Next another situation is considered. If one J2 is ex- 
changed by Ji in the alternating chain, the system has a 
configuration 



• J\ J% J\ Ji J\ 02 J\ 



(17) 



We define A& = for this configuration. Shifting a posi- 
tion of a J\ by one, the system has two J\J\ pairs and 
has a configuration, 



• • • J\J\JiJ\J\ 



(18) 



We define this distance between J\ J\ and J\ J\ as Ab = 1 . 
Furthermore, Ab = 2 is defined for the configuration jl9| ) 
and etc. 



■ J1J1J2J1J2J1J1 



(19) 



The ground state energies for a PBC chain of L = 24 are 
shown as a function of A b in Fig. ^ (b). We also study 
the interaction between J2J2 and J2J2, where we define 
the distance between the impurities in the same way. We 
again found that the ground state energy increases as the 
distance becomes larger. Thus it has been found that an 
attractive force acts between the impurities regardless of 
their type. 

If we allow positions of impurities to move, the distance 
between impurities would be distributed in the canoni- 
cal distribution according to the interaction between the 
impurities in the thermal equilibrium at a given temper- 
ature as has demonstrated in the previous paper p7[ . 



V. RESPONCE OF THE LOCAL MAGNETIC 
STRUCTURE TO THE DYNAMICAL FIELD 

In this section we investigate dynamical response of the 
local magnetic structure (effective S — 1/2 spin) induced 
by a bond impurity to a time-dependent magnetic field 
with a random noise in the dynamical model (^]). 

First let us consider the response of the magnetization 
of a free S = 1/2 spin to a sweeping field. The Hamilto- 
nian is given by 



H{t) = 2TS x -(-H + ct)S z . 



(20) 



This system is two-level system whose energy levels are 
shown in Fig. |as a function of H , where the avoided 
level crossing occurs near H = 0. If Hq 3> T > 0, the 
ground state consists primarily of the M z = — 1/2 state 
at t = 0. The probability for the system to end up in 
the M z = 1/2 state (i.e., the probability to change its 
magnetization) at t — 00 is given by Eq. ( |l0| ) 

Next we consider the case of a local magnetic structure 
(effective S = 1/2) on a lattice. When the noise {hi(t)} 
does not exist, dynamics of the system is generally the 
same as that of a free S = 1/2 spin as far as we concern 
the states in an S = 1/2 space. It is easily understood as 
follows. Adopting { |+), |— )} as the basis set, the matrix 
representation of Eq. (EOh is 



n 



-f r 
r f 



(21) 



Noting that ]T\ J M+1 S< ■ S l+1 and 2r ^ Sf - H(t) £\ Sf 
commute, the matrix representation of Eq. (g) in the 
total S =1/2 space is given by 



n 



H_ 

2 



const. 



r 



4f + const. 



(22) 



adopting {|5 tot = 1/2, M z = 1/2), IS" * = 1/2, M z = 
— 1/2)} as the basis set. The constant in Eq. ( |22| ) does 
not depend on H(t). Thus the dynamics is independent 
of the number of sites and the combination of { Jij}- 

In the experimental situation, however, noise usually 
has an influence on the system. If a random noise is 
applied to each site individually, the term — ^2^(11 (t) + 
hi(t))S? does not commute with (S tot ) 2 = (£ s *) 2 ; and 
therefore the dynamics would be changed when the in- 
teraction {Jij} and the system size are varied. 

Here for the study of dynamical properties we adopt 
the minimum model of the type (d) because we can treat 
a limited number of spins in the scheme of Eq. @ . This 
model consists of five sites and we take J\ = 1.0 and 
J2 = 0.5. This is the minimum bond-alternating model 
with a J2J2 bond impurity. This system is looked on as 
an effective S =1/2 spin. 

It would be expected that the effect of noise is reduced 
as the system size becomes large and a local magnetic 
structure (effective S = 1/2 spin) consisting of several 
number of spins is less sensitive to such a random noise 
than a free 5 = 1/2 spin. In order to study effects of the 
noise we investigate the following properties. 

First we study the broadening of the level due to the 
noise. A energy gap and a sweeping rate at the level 
crossing point are important factors to determine the 
transition probability under a sweeping uniform field (see 
Eq. (p^|)). It would be useful to investigate how the noise 
influences the energy gap at the level crossing point. For 
a single spin the energy gap is given by 



ae = V(2r) 2 + {h{t)y 



(23) 



for H(t) = 0. Thus using the distribution of h(t) (a 
gaussian distribution with the variance A 2 /2^f in Eq. 
(||)), we can calculate the mean (AE) and the width 
SE = yf((AE) 2 ) - (AE) 2 , which are listed in Table | 
For the local magnetic structure, we calculate the energy 
gap in a noise by a perturbation method. That gives 



AE = 2r 



W'IW)| a 



E. 



(0) 



E 



(0) 



(24) 



+ contributions from higher levels, 



where Ei and (j)^ 1 (E 2 V> and 



(0) / 77,(0) 



are the eigenstate 



of the ground state (the first exited state) of the non 
perturbed hamiltonian. The term V\ is hiSf . Noting 
that 
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(K^ 0) l^l4 0) )l 2 ) = E^ 2 )K^ 0) l^l^ 0) )l 2 , (25) 

i 

we can calculate (AE) from the distribution of {hi}, 
which are also listed in Table | for L = 5 and L = 9. 
For large values of A, we obtain the energies by diago- 
nalizing the hamiltonian Eq. (Q) and obtained (AE) and 
SE numerically from 500 samples of {hi}, which are also 
listed in Table |. Thus we find the noise causes almost 
the same effect on the broadening of the energy levels at 
H(t) — in both systems, i.e., a single spin and local 
magnetic structures, which is not in accordance with the 
above expectation. 

Second, we investigate the dynamical properties of the 
both systems. In particular, we study the time evolu- 
tion of magnetization under the sweeping field. We set 
parameters V = 0.02, H = 0.5, c = 0.0005, i max =2000, 
and dt = 0.01. In this parameter set the probability p in 
Eq. ( pp| ) is nearly 1. As was mentioned, the responses 
of the magnetization are the same in both systems when 
a random noise is not applied. Time evolution without 
noise is shown by thin dotted lines in Figs. || (a)-(c). 

Next we investigate the dynamics with the random 
field (Eq. (H)). We observe the two cases changing the 
amplitude A of noise (Eq. (0)); (a) A = 0.01 and (b) 
A = 0.02. The value of 7 is fixed at 0.1. The lowest two 
energy levels of the system with a random noise {hi(t)} 
for case (b) are illustrated as a function of time in Fig. 

Figure H (a) shows M(t) in the noise of a small am- 
plitude (A = 0.01). The transition probabilities for a 
local magnetic structure (effective 5 = 1/2 spin) and an 
5 = 1/2 free spin are very close to each other at this 
small amplitude and the probabilities are reduced by a 
little amount from the probability of the pure system 
(i.e., hi(t) = 0). Increasing the amplitude to A = 0.02 
(Fig. H (b)), reduction of M(t)s increases. Furthermore 
we find that M(t) for the local magnetic structure re- 
mains larger than that of a free spin. The errorbar shows 
the first standard deviation of the distribution of {M(t)}. 
Here we find that the distribution is rather wide but the 
mean of the distribution is definitely different. The re- 
duction of M(t) by field fluctuation corresponds to the 
reduction of p in Eq. (|To|). The reduction of p has al- 
ready pointed out for a single spin by Y. Kayanuma and 
H. Nakayama ^9|. Our results are qualitatively consis- 
tent with their results. Furthermore, the present obser- 
vation indicates that a local magnetic structure induced 
by an impurity is less sensitive to the noise. This feature 
may be due to this bulky structure of the local magnetic 
structure. From this observation we may expect that lo- 
cal magnetizations are easier to manipulate than single 
spins in a field with noise. Detail analysis for the noise 
dependence will be reported elsewhere. 



VI. SUMMARY AND DISCUSSION 

In the bond-alternating (5 = 1/2) Heisenberg anti- 
ferromagnetic chain, we studied properties of locally in- 
duced magnetic structures by inhomogeneities of the lat- 
tice such as the defect of alternation or the edge of the 
lattice. Because of the gapful nature, the role of the inho- 
mogeneities is similar to that in the Haldane systems and 
is very different from that in the uniform chain. In the 
uniform chain we found that the bond impurities divide 
the system into domains. We showed that a local mag- 
netic structure induced by a bond impurity can be looked 
on as an effective 5 = 1/2 spin. The interaction between 
the local magnetic moments decays exponentially in the 
present model and the local magnetic structures behave 
almost independently even at very low temperatures. 

We also studied the force acting between these defects 
of alternation. There are several configurations of de- 
fects. We investigated the forces between J\J\ and J\J\, 
between J 2 Ji and JiJ2-, and between J\J\ and JiJ2- Fur- 
thermore, there are two ty pes of separations of defects, 
i.e., the type of ( |l5| ) and (|l7|). It turned out that the 
forces are attractive for all cases. 

We also studied the dynamical response of the magne- 
tization to a sweeping field. The induced local moment 
behaves as a single spin even in dynamical property as 
far as the uniform field is applied. We considered what 
property is different between the local magnetic struc- 
ture (effective 5=1/2 spin) induced by an impurity and 
a free 5=1/2 spin. In realistic situation there exists 
noise which is not necessarily uniform but different from 
site to site. It has been found that even when noise is 
applied individually at each site, the total effect of the 
noise on the broadening of the energy levels is almost the 
same as that of the single spin. However, the dynamical 
response of the magnetization to a time-dependent mag- 
netic field with a random noise is found to be different 
from that of a single spin. That is, the local magnetic 
structure is found to be more robust against noise applied 
at each site individually. 

The present study is the first attempt treating the dy- 
namics of such a local magnetic structure. We hope that 
this work provides a basic information for the manipula- 
tion of microscopic or nanoscale magnetic devices in the 
future. We also hope that the present study helps to an- 
alyze the magnetic properties at low temperatures such 
as observed by NMR measurement. 
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APPENDIX A: THE CORRELATION LENGTH 
OF BOND- ALTERNATING CHAIN (S = 1/2) 

We determined the correlation length of HAF bond- 
alternating chain (S — 1/2) in the ground state as a 
function of the ratio Jx/J^- We calculated the correla- 
tion function for various ratios (J\/ J2) and estimated the 
correlation length using the relation 



ln(S? S?) = constant - 



(Al) 



where denotes the distance between site i and site j, 
and £ denotes the correlation length of the system. We 
treated periodic chains with 60 sites and simulations were 
performed at T — 0.01 in the M z = space. We show 
the correlation lengths for various ratios (J1/J2) in Table 
y. Here we adopt the length of a pair of bonds Ji and J2 
to be a unit length, i.e., this unit is double of the bond 
length. 
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FIG. 1. Magnetization profiles {ra;} of the models (a)-(d) 
at T — 0.01 in the M z = 1/2 space, (a) and (c) contain 63 
sites (L — 63) and (b) and (d) contain 65 sites (L — 65) . 
The diamonds denote the strength of bonds { J^i+i}, those at 
high positions denote Ji = 1.3 and those at the low positions 
denote J2 = 0.7. Details are shown in text. 



FIG. 2. Summation of magnetization per site from the 
left edge site for the model (b). 



FIG. 3. Summation of magnetization per site from the left 
edge site for a bond-alternating chain of the type (b) in Fig. 
|l] with L = 21. The strong and weak bonds are Ji = 2 and 
J2 = 0.5, respectively. 



FIG. 4. Distribution of M z for the model (b) in Fig. [j]. 



FIG. 5. Magnetization profile {m^} for a model which has 
five positions of induced moments in the M z = 1/2 space at 
T — 0.01 with L — 101. The diamonds denote the strength of 
bonds {J^i+i}, those at high positions denote Ji — 1.3 and 
those at the low positions denote J2 = 0.7. 
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FIG. 6. Ground state energy as a function of (a) A a and 
(b) A;, obtained by an exact diagonalization [L — 24, PBC). 



FIG. 7. Avoided level crossing. 

FIG. 8. Time evolution of magnetization. Thin dotted 
line denotes the magnetization without noise. Thick line de- 
notes the magnetization of a local magnetic structure (effec- 
tive S — 1/2 spin) of L = 5 with a random noise (averaged 
over 500 samples {hi}). Thick dotted line denotes the mag- 
netization of a free S — 1/2 spin with a random noise (av- 
eraged over 500 samples {hi}), (a) A = 0.01,7 = 0.1, (b) 
A = 0.02,7 = 0.1. 

FIG. 9. Energy structure of the two level system with a 
random noise for A = 0.02, 7 = 0.1. 



TABLE I. Average energy gaps (AE) and their widths SE 
for several values of A/y/2j. (1), (5), and (9) mean a single 
free spin, the minimum model consisting of 5 spins, and the 
next minimum model consisting of 9 spins of the type (d) 
in Fig. [jj, respectively. The subscript p denotes that the 
quantities are obtained by the perturbation method and the 
subscript d denotes that the quantities are obtained by the 
exact diagonalization method. 



TABLE II. Correlation length £ as a function of the ratio 
of the J1/J2 in the bond-alternating chain (S = 1/2) 
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Table I 





AE d (l) 


SE d (l) 


A£ p (5) 


A£ d (5) 


5E d {5) 


A£ P (9) 


AE d {9) 


6E d (9) 


0.01 


0.0412 


0.0016 


0.0413 


0.0412 


0.0016 


0.0413 


0.0413 


0.0017 


0.02 


0.0445 


0.0056 


0.0450 


0.0443 


0.0054 


0.0450 


0.0446 


0.0058 


0.03 


0.0486 


0.0105 


0.0513 


0.0492 


0.0111 


0.0514 


0.0492 


0.0112 


0.04 


0.0542 


0.0157 


0.0601 


0.0544 


0.0169 


0.0602 


0.0546 


0.0171 


0.05 


0.0609 


0.0222 


0.0714 


0.0591 


0.0213 


0.0716 


0.0605 


0.0233 



Table II 



Jl/J2 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


1.9 


2.0 




2.02 


1.51 


1.31 


1.18 


1.01 


0.93 


0.84 


0.80 


0.77 



